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Abstract 

A search for the decays B® — > /i + jjT and B° — >■ \x~ is performed with 0.37 fb _1 of 
collisions at a/s = 7 TeV collected by the LHCb experiment in 2011. The upper 
limits on the branching fractions are B(B° S — > y^jJT) < 1.6 x 10~ 8 and B(B° — > 

< 3.6 x 10~ 9 at 95 % confidence level. A combination of these results with the LHCb limits 
obtained with the 2010 dataset leads to B(B° S < 1.4 x 10" 8 and B(B° 

< 3.2 x 10~ 9 at 95% confidence level. 
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1 Introduction 



Measurements of low-energy processes can provide indirect constraints on particles that 
are too heavy to be produced directly. This is particularly true for Flavour Changing 
Neutral Current (FCNC) processes which are highly suppressed in the Standard Model 
(SM) and can only occur through higher-order diagrams. The SM predictions for the 
branching fractions of the FCNC decaya^j-B 9 — > and B° — > are B(B® — > p + p~) 

= (3.2±0.2) x 10" 9 and B(B° = (0.10±0.01) x 1(T 9 |1|. However, contributions 

from new processes or new heavy particles can significantly enhance these values. For 
example, within Minimal Supersymmetric extensions of the SM (MSSM), in the large 
tan/3 regime, B(B® — > ^ + fi~) is found to be approximately proportional to tan 6 /3 [2], 
where tan is the ratio of the vacuum expectation values of the two neutral CP-even 
Higgs fields. The branching fractions could therefore be enhanced by orders of magnitude 
for large values of tan f3. 

The best published limits from the Tevatron are B(B® — > < 5.1 x 1CT 8 

at 95% confidence level (CL) by the DO collaboration using 6.1 fb _1 of data [3], and 
B(B° -» fM+fj,-) < 6.0 x 10~ 9 at 95% CL by the CDF collaboration using 6.9 fb^ 1 of 
data pi. In the same dataset the CDF collaboration observes an excess of B° s — > [i + 
candidates compatible with B(B® — > = (1.8i 'g) x 10~ 8 and with an upper limit of 

B(B° S ->- < 4.0 x 10" 8 at 95% CL. The CMS collaboration has recently published 

B(B° S ->• fi+fi') < 1.9 x 10- 8 at 95% CL and B(B° < 4 -6 x 10" 9 at 95% CL using 

1.14 fb _1 of data pi. The LHCb collaboration has published the limits [H] B(B® — > 
< 5.4 x 10^ 8 and~B(B° ^ + ^~) < 1.5 x 10^ 8 at 95% CL based on about 37 pb" 1 of 
integrated luminosity collected in the 2010 run. 

This Letter presents an analysis of the data recorded by LHCb in the first half of 2011 
which correspond to an integrated luminosity of ~ 0.37 fb -1 . The results of this analysis 
are then combined with those published from the 2010 dataset. 

2 The LHCb detector 

The LHCb detector [7] is a single-arm forward spectrometer designed to study production 
and decays of hadrons containing b or c quarks. The detector consists of a vertex loca- 
tor (VELO) providing precise locations of primary pp interaction vertices and detached 
vertices of long lived hadrons. 

The momenta of charged particles are determined using information from the VELO 
together with the rest of the tracking system, composed of a large area silicon tracker 
located before a warm dipole magnet with a bending power of ~ 4 Tm, and a combination 
of silicon strip detectors and straw drift chambers located after the magnet. Two Ring 
Imaging Cherenkov (RICH) detectors are used for charged hadron identification in the 
momentum range 2-100 GeV/c. Photon, electron and hadron candidates are identified by 
electromagnetic and hadronic calorimeters. A muon system of alternating layers of iron 

inclusion of charged conjugated processes is implied throughout. 
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and drift chambers provides muon identification. The two calorimeters and the muon 
system provide the energy and momentum information to implement a first level (LO) 
hardware trigger. An additional trigger level (HLT) is software based, and its algorithms 
are tuned to the experimental operating condition. 

Events with a muon final states are triggered using two LO trigger decisions: the 
single-muon decision, which requires one muon candidate with a transverse momentum 
Pr larger than 1.5 GeV/c, and the di-muon decision, which requires two muon candidates 
with transverse momenta pp,i and px,2 satisfying the relation a/Pt,i • Pt,2 > 1-3 GeV/c. 

The single muon trigger decision in the second trigger level (HLT) includes a cut on 
the impact parameter (IP) with respect to the primary vertex, which allows for a lower 
Pt requirement (px > 1-0 GeV/c, IP > 0.1 mm). The di-muon trigger decision requires 
muon pairs of opposite charge with > 500 MeV/c, forming a common vertex and 
with an invariant mass > 4.7 GeV/c 2 . A second trigger decision, primarily to select 
J/if> events, requires 2.97 < m w < 3.21 GeV/c 2 . The remaining region of the di-muon 
invariant mass range is also covered by trigger decisions that in addition require the 
di-muon secondary vertex to be well separated from the primary vertex. 

Events with purely hadronic final states are triggered by the LO trigger if there is a 
calorimeter cluster with transverse energy > 3.6 GeV. Other HLT trigger decisions 
select generic displaced vertices, providing high efficiency for purely hadronic decays. 



3 Analysis strategy 

Assuming the branching fractions predicted by the SM, and using the bb cross-section 
measured by LHCb in the pseudorapidity interval 2 < rj < 6 and integrated over all 
transverse momenta of o b i = 75 ± 14/xb j8], approximately 3.9 B° — > /i + [T and 0.4 
B° events are expected to be triggered, reconstructed and selected in the analysed 

sample embedded in a large background. 

The general structure of the analysis is based upon the one described in Ref. [6]. First 
a very efficient selection removes the biggest amount of background while keeping most of 
the signal within the LHCb acceptance. The number of observed events is compared to 
the number of expected signal and background events in bins of two independent variables, 
the invariant mass and the output of a multi-variate discriminant. The discriminant is 
a Boosted Decision Tree (BDT) constructed using the TMVA package M. It supersedes 
the Geometrical Likelihood (GL) used in the previous analysis |6) as it has been found 
more performant in discriminating between signal and background events in simulated 
samples. No data were used in the choice of the multivariate discriminant in order not to 
bias the result. 

The combination of variables entering the BDT discriminant is optimized using sim- 
ulated events. The probability for a signal or background event to have a given value of 
the BDT output is obtained from data using B9, — > h + h ~ candidates (where ' can be 
a pion or a kaon) as signal and sideband B? s ^ — > fi + candidates as background. 



The invariant mass line shape of the signals is described by a Crystal Ball function 10 
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with parameters extracted from data control samples. The central values of the masses are 
obtained from B° — > K + tc~ and B® — > K + K~ samples. The B° s and B° mass resolutions 
are estimated by interpolating those obtained with di-muon resonances ( J/i/j, ?p(2S) and 
T(1S,2S,3S)) and cross-checked with a fit to the invariant mass distributions of both 
inclusive B9^ — » h + h ~ decays and exclusive B° — > K + tt~ decays. The central values of 
the masses and the mass resolution are used to define the signal regions. 

The number of expected signal events, for a given branching fraction hypothesis, is 
obtained by normalizing to channels of known branching fractions: B + — > J/ipK + , B®— > 
J/ip4> and B° — > K + n~ . These channels are selected in a way as similar as possible to 
the signals in order to minimize the systematic uncertainty related to the different phase 
space accessible to each final state. 

The BDT output and invariant mass distributions for combinatorial background events 
in the signal regions are obtained using fits of the mass distribution of events in the mass 
sidebands in bins of the BDT output. 

The two-dimensional space formed by the invariant mass and the BDT output is 
binned. For each bin we count the number of candidates observed in the data, and 
compute the expected number of signal events and the expected number of background 
events. The binning is unchanged with respect to the 2010 analysis [6]. The compatibility 
of the observed distribution of events in all bins with the distribution expected for a given 
branching fraction hypothesis is computed using the CL S method jll], which allows a 
given hypothesis to be excluded at a given confidence level. 



4 Selection 

The B9 g s — > fi + fi~ selections require two muon candidates of opposite charge. Tracks are 
required to be of good quality and to be displaced with respect to any primary vertex. 
The secondary vertex is required to be well fitted (x 2 /nDoF < 9) and must be sepa- 
rated from the primary vertex in the forward direction by a distance of flight significance 
(L/a(L)) greater than 15. When more than one primary vertex is reconstructed, the one 
that gives the minimum impact parameter significance for the candidate is chosen. The 
reconstructed candidate has to point to this primary vertex (IP/cr(IP) < 5). 

Improvements have been made to the selection developed for 2010 data [6]. The RICH 
is used to identify kaons in the B® — > J/i(j<j) normalization channel and the Kullback- 



Leibler (KL) distance 12 is used to suppress duplicated tracks created by the recon- 
struction. This procedure compares the parameters and correlation matrices of the recon- 
structed tracks and where two are found to be similar, in this case with a symmetrized 
KL divergence less than 5000, only the one with the higher track fit quality is considered. 

The inclusive B9s — > h + h~ sample is the main control sample for the determina- 
tion from data of the probability distribution function (PDF) of the BDT output. This 
sample is selected in exactly the same way as the B9, — > fi + fi~ signals apart from the 
muon identification requirement. The same selection is also applied to the B° — > K + n~ 
normalization channel. 
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The muon identification efficiency is uniform within ~ 1% in the considered phase 
space therefore no correction is added to the BDT PDF extracted from the B^ — > h + h~ 
sample. The remaining phase space dependence of the muon identification efficiency 
is instead taken into account in the computation of the normalization factor when the 
B° — > K + tt~ channel is considered. 

The J/ip — > fifi decay in the B + — > J/ipK + and B® — > J/ifxft normalization channels 
is selected in a very similar way to the B9\ — > channels, apart from the pointing 

requirement. K ± candidates are required to be identified by the RICH detector and to 
pass track quality and impact parameter cuts. 

To avoid pathological events, all tracks from selected candidates are required to have 
a momentum less than ITeV/c. Only B candidates with decay times less than 5r B o , 

( s ) 



where t b o is the B lifetime 13 , are accepted for further analysis. Di-muon candidates 



coming from elastic di-photon production are removed by requiring a minimum transverse 
momentum of the B candidate of 500MeV/c. 



5 Determination of the mass and BDT distributions 

The variables entering the BDT discriminant are the six variables used as input to the 
GL in the 2010 analysis plus three new variables. The six variables used in the 2010 anal- 
ysis are the B lifetime, impact parameter, transverse momentum, the minimum impact 
parameter significance (IP/a(IP)) of the muons, the distance of closest approach between 
the two muons and the isolation of the two muons with respect to any other track in the 
event. The three new variables are: 



1. the minimum of the two muons; 



2. the cosine of the angle between the muon momentum in the B rest frame and the 
vector perpendicular to the B momentum and the beam axis: 

Pt,b («fy M /2) 

where /ii labels one of the muons and is the reconstructed B candidate mass^j 

3. the B isolation 14 

where Pt{B) is the B transverse momentum with respect to the beam line 
and the sum is over all the tracks, excluding the muon candidates, that satisfy 



3 As the B is a (pseudo)-scalar particle, this variable is uniformely distributed for signal candidates 
while is peaked at zero for bb — > /i + X background candidates. In fact, muons from semi-leptonic decays 
are mostly emitted in the direction of the 6's and, therefore, lie in a plane formed by the B momentum 
and the beam axis. 
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a/ 5r] 2 + 5(j) 2 < 1.0, where 5rj and 50 denote respectively the difference in pseudo- 
rapidity and azimuthal angle between the track and the B candidate. 

The BDT output is found to be independent of the invariant mass for both signal and 
background and is defined such that the signal is uniformly distributed between zero and 
one and the background peaks at zero. The BDT range is then divided in four bins of 
equal width. The BDT is trained using simulated samples {B^ fi + for signals and 
bb — > fi + fi~X for background where X is any other set of particles) and the PDF obtained 
from data as explained below. 

5.1 Combinatorial background PDFs 

The BDT and invariant mass shapes for the combinatorial background inside the signal 
regions are determined from data by interpolating the number of expected events using 
the invariant mass sidebands for each BDT bin. The boundaries of the signal regions are 
defined as m^o ± 60 MeV/c 2 and mgo ±60 MeV/c 2 and the mass sidebands as [m^o — 
600 MeV/c 2 , m B o - 60 MeV/c 2 ] and [m B o + 60 MeV/c 2 , m B o + 600 MeV/c 2 ]. 

Figure [T] shows the invariant mass distribution for events that lie in each BDT output 
bin. In each case the fit model used to estimate the expected number of combinatorial 
background events in the signal regions is superimposed. 

Aside from combinatorial background, the low-mass sideband is potentially polluted 
by two other contributions: cascading b — > cfiu — > fi/iX decays below 4900MeV/c 2 and 
peaking background from B? a ■> — ^ h + h~ candidates with the two hadrons misidentified 
as muons above 5000 MeV/c 2 . To avoid these contaminations, the number of expected 
combinatorial background events is obtained by fitting a single exponential function to 
the events in the reduced low-mass sideband [4900, 5000] MeV/c 2 and in the full high- 
mass sideband. As a cross-check, two other models, a single exponential function and the 
sum of two exponential functions, have been used to fit the events in different ranges of 
sidebands providing consistent background estimates inside the signal regions. 

5.2 Peaking background PDFs 

The peaking backgrounds due to B9, — > h + h'~ events in which both hadrons are misiden- 
tified as muons have been evaluated from data and simulated events to be N B o = 1.0 ±0.4 
events and N B o = 5.0 ± 0.9 events within the two mass windows and in the whole BDT 
output range. The mass line shape of the peaking background is obtained from a simu- 
lated sample of doubly-misidentified B?* — > h + h ~ events and normalized to the number 
of events expected in the two search windows from data, N B o and N B o. The BDT PDF 
of the peaking background is assumed to be the same as for the signal. 

5.3 Signal PDFs 

The BDT PDF for signal events is determined using an inclusive B^ — > h + h'~ sam- 
ple. Only events which are triggered independently on the signal candidates have been 
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Figure 1: Distribution of the invariant mass for events in each BDT output bin. 

The curve shows the model used to fit the sidebands and extract the expected number 
of combinatorial background events in the B® and B° signal regions, delimited by the 
vertical dotted orange and dashed green lines respectively. Only events in the region in 
which the line is solid have been considered in the fit. 



considered (TIS events). 

The number of B^ — > h + h~ signal events in each BDT output bin is determined by 



fitting the hh! invariant mass distribution under the mass hypothesis 15 . Figure \i 



shows the fit to the mass distribution of the full sample and for the three highest BDT 
output bins for B9 , — > h + h ~ TIS events. The B9 g , — > h + h ~ exclusive decays, the 
combinatorial background and the physical background components are drawn under the 
fit to the data; the physical background is due to the partial reconstruction of three-body 
B meson decays. 

In order to cross-check this result, two other fits have been performed on the same 
dataset. The signal line shape is parametrized either by a single or a double Crys- 



tal Ball function 10 , the combinatorial background by an exponential function and 
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In addition, exclusive B9s — > 



the physical background by an ARGUS function 
ir~K + , 7r~7r + , K~K + channels, selected using the K — ir separation capability of the 
RICH system, are used to cross-check the calibration of the BDT output both using 
the tt~K + , 7r~7r + , K~K + inclusive yields without separating B and B° s and using the 
B° — > K + 7i~ exclusive channel alone. The maximum spread in the fractional yield ob- 
tained among the different models has been used as a systematic uncertainty in the signal 
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m (MeV/c 2 ) <n m (MeV/c 2 ) 



Figure 2: Invariant mass distributions of B9x — > h + h'~ candidates in the mass 
hypothesis for the whole sample (top left) and for the samples in the three highest bins 
of the BDT output (top right, bottom left, bottom right). The B^ — > h + h ~ exclusive 
decays, the combinatorial background and the physical background components are drawn 
under the fit to the data (solid blue line). 



BDT PDF. The BDT PDFs for signals and combinatorial background are shown in Fig. |3j 
The invariant mass shape for the signal is parametrized as a Crystal Ball function. The 
mean value is determined using the B° — > K + n~ and B° s — > K + K~ exclusive channels 
and the transition point of the radiative tail is obtained from simulated events p\. The 
central values are 

m B o = 5358.0 ± 1.0 MeV/c 2 , 
m B o = 5272.0 ± 1.0 MeV/c 2 . 

The measured values of m B o and m B o are 7 — 8 MeV/c 2 below the PDG values 13 
due to the fact that the momentum scale is uncalibrated in the dataset used in this 
analysis. The mass resolutions are extracted from data with a linear interpolation between 
the measured resolution of charmonium and bottomonium resonances decaying into two 
muons: J/ip, ip(2S), T(1S), T(2S) and T(35). The mass line shapes for quarkonium 
resonances are shown in Fig. |4j Each resonance is fitted with two Crystal Ball functions 
with common mean value and common resolution but different parameterization of the 
tails. The background is fitted with an exponential function. 

The results of the interpolation at the m B o and m B o masses are 

a(m B o) = 24.6±0.2 (stat) ±1.0 (syst) MeV/c 2 , 
a(m B o) = 24.3±0.2 (stat) ±1.0 (syst) MeV/c 2 . 
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Figure 3: BDT probability distribution functions of signal events (solid squares) and 
combinatorial background (open circles): the PDF for the signal is obtained from the 
inclusive sample of TIS — > h + h~ events, the PDF for the combinatorial background 
is obtained from the events in the mass sidebands. 
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Figure 
(9-11) 



4: Di-muon invariant mass spectrum in the ranges (2.9 - 3.9) GeV/c 2 (left) and 
MeV/c 2 (right). 



This result has been checked using both the fits to the B9 \ — > h + fi inclusive decay line 
shape and the B° — > K + n~ exclusive decay. The results are in agreement within the 
uncertainties. 



S 



6 Normalization 



To estimate the signal branching fraction, the number of observed signal events is 
normalized to the number of events of a channel with a well known branching frac- 
tion. Three complementary normalization channels are used: B + — > J/ip({i + fi~)K + , 
B° s J/0 fi~)<f)(K + K~) and B° — » K + ti~ . The first two channels have similar trigger 
and muon identification efficiencies to the signal but different number of particles in the 
final state. The third channel has a similar topology but is selected by different trigger 
lines. 

The numbers of B® — > fi + fi~ and B° — > candidates are translated into a branch- 

ing fractions (£>) using the equation 



,REC SEL|REC TRIG|SEL , N 

K? _ Y3 s , fc norm fc norm tnorm ./norm ^ D (s)^^ V _ ^ norm AT /o\ 

° rm rec SEL|REC TRIG|SEL f X jj^ ~ Q B^^- X Bf> B) ^+ , W 



-sig 



where /^( s ) and / norm are the probabilities that a 6 quark fragments into a TiPs and 
into the b hadron involved for the chosen normalization mode. LHCb has measured 



fs/fd — 0-267^q;q2o 17 • -^norm is the branching fraction and iV norm is the number of se- 
lected events of the normalization channel. The efficiency is the product of three factors: 
£ rec j g ^ e reconstruction efficiency of all the final state particles of the decay including the 
geometric acceptance of the detector; e SEL l REC i s the selection efficiency for reconstructed 
events; g TRIG l SEL is the trigger efficiency for reconstructed and selected events. The sub- 
script (sig, norm) indicates whether the efficiency refers to the signal or the normalization 
channel. Finally, a^°o rm _ s>M + At - is the normalization factor (or single event sensitivity) and 

Nro .,,+ ,,- the number of observed signal events. 

For each normalization channel N noxm is obtained from a fit to the invariant mass 
distribution. The invariant mass distributions for reconstructed B + — > J / tpK + and B® — > 
J/ipcj) candidates are shown in Fig. |5j while the B° — > K + ir~ yield is obtained from the 
full B9 g -, " fit as shown in the top left of Fig. 2 

The numbers used to calculate the normalization iactors are summarized in Table [TJ 
A weighted average of the three normalization channels, assuming the tracking and trig- 
ger uncertainties to be correlated between the two J/ip normalization channels and the 
uncertainty on fd/ f s to be correlated between the B + — > J/ipK + and B° — > K + n~, gives 



«^ +M - = (8.38 ± 0.74) xlO- lu , 
<°"W = (2-20 ± 0.11) xlO- 10 . 
These normalization factors are used to determine the limits. 



7 Results 

The results for B® — > and B° — > fi + fi~ are summarized in Table [2] and Table [3 

respectively and in each of the bins the expected number of combinatorial background, 
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Table 1: Summary of the quantities and their uncertainties required to calculate the 



normalization factors 



'MA* 



for the three normalization channels considered. The 



branching fractions are taken from Refs. 13 18 . The trigger efficiency and the number 



of B° — > K + n candidates correspond to TIS events. 



REC SEL | REC TRIG|SEL 

K £ norm e norm £ norm AT IlOrm _,norm 

° eRB C e SEL|REC £ TRIG|SEL JV norm Ut go _^ ^ + ^_ a B°^+p- 

sig sig sig 

(xio- 5 ) (xitr 10 ) (xio- 9 ) 

B+^J/ipK + 6.01±0.21 0.48±0.014 0.95±0.01 124518±2025 2.23±0.11 0.83±0.08 

B s ° -> J/^ 3.4±0.9 0.24±0.014 0.95±0.01 6940±93 2.96±0.84 1.11±0.30 

B" K+n- 1.94±0.06 0.86±0.02 0.049±0.004 4146±608 1.98±0.34 0.74±0.14 



peaking background, signal events, with the SM prediction assumed, is shown together 
with the observations on the data. The uncertainties in the signal and background PDFs 
and normalization factors are used to compute the uncertainties on the background and 
signal predictions. 

The two dimensional (mass, BDT) distribution of selected events can be seen in Fig. [6] 
The distribution of the invariant mass in the four BDT bins is shown in Fig. [7] for B® — >■ 
and in Fig. [8] for B° — > fi + fi~ selected candidates. 

The compatibility of the distribution of events inside the search window in the invari- 
ant mass-BDT plane with a given branching fraction hypothesis is evaluated using the 
CL S method |11|. This method provides three estimators: CL s+b , a measure of the com- 
patibility of the observed distribution with the signal and background hypotheses, CLt,, a 
measure of the compatibility with the background-only hypothesis and CL S , a measure of 
the compatibility of the observed distribution with the signal and background hypotheses 
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Figure 6: Distribution of selected di-muon events in the invariant mass-BDT plane. The 
orange short-dashed (green long-dashed) lines indicate the ±60 MeV/ c 2 search window 
around the mean B° s (B°) mass. 

normalized to the background-only hypothesis. 

The expected CL S values are shown in Fig. [9] for — >■ and for B° — > 

as dashed black lines under the hypothesis that background and SM events are observed. 
The shaded areas cover the region of ±lcr of compatible observations. The observed values 
of CL S as a function of the assumed branching ratio is shown as dotted blue lines on both 
plots. 

The expected limits and the measured limits for B® — > fi + fi~ and B° — > fi + fi~ at 90 % 
and 95 % CL are shown in Table [4] and Table [HJ respectively. For the B® — > decay, 
the expected limits are computed allowing the presence of B° s — > ji + fi~ events according 
to the SM branching fraction. For the B° — > fi + fi~ decay the expected limit is computed 
in the background-only hypothesis and also allowing the presence of B° — > \x + \x~ events 
with the SM rate: the two results are identical. In the determination of the limits, the 
cross-feed of B® — > n + fi~ (B° — > fi + n~) events in the B° (B®) mass window has been 
taken into account assuming the SM rates. 

The observed CLb values are shown in the same tables. The comparison of the observed 
distribution of events with the expected background distribution results in a p-value 
(1 - CL b ) of 5% for the B° s ->• and 32% for the B° decay. For the 

B° s — > ji + \i~ decay, the probability that the observed events are compatible with the sum 
of expected background events and signal events according to the SM rate is measured 
by 1— CL s+b and it is 33%. 

The result obtained in 2011 with 0.37 fb _1 has been combined with the published 
result based on ~ 37 pb" 1 |6|. The expected and observed limits for 90 % and 95 % CL 
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Figure 7: Distribution of selected di-muon events in the B® — » mass window for 

the four BDT output bins. The black dots are data, the light grey histogram shows 
the contribution of the combinatorial background, the black filled histogram shows the 
contribution of the B9 x — > h + h~ background and the dark grey filled histogram the 
contribution of B° s — > fi + n~ signal events according to the SM rate. The hatched area 
depicts the uncertainty on the sum of the expected contributions. 



for the combined results are shown in Table [4] for the B° s — > fi + fi decay and in Table [5] 
for the B° — > fi + fi~ decay. 



8 Conclusions 

With 0.37 ftT 1 of integrated luminosity, a search for the rare decays B® — > fi + and 
B° — > has been performed and sensitivities better than the existing limits have been 

obtained. The observed events in the B° s and in the B° mass windows are compatible 
with the background expectations at 5% and 32% confidence level, respectively. For the 
B° s fi + fi~ decay, the probability that the observed events are compatible with the sum 
of expected background events and signal events according to the SM rate is 33%. The 
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Figure 8: Distribution of selected di-muon events in the B° — > /i + yT mass window for 
the four BDT output bins. The black dots are data, the light grey histogram shows 
the contribution of the combinatorial background, the black filled histogram shows the 
contribution of the B9\ — )■ h + h~ background and the dark grey filled histogram shows 
the cross-feed of B® — > fi + fi~ events in the B° mass window assuming the the SM rate. 
The hatched area depicts the uncertainty on the sum of the expected contributions. 
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Figure 9: CL S as a function of the assumed B. Expected (observed) values are shown by 
dashed black (dotted blue) lines. The expected CL S values have been computed assuming 
a signal yield corresponding to the SM branching fractions. The green (grey) shaded areas 
cover the region of ±lcr of compatible observations. The measured upper limits at 90% 
and 95% CL are also shown. Left: B® — > ji + jj," , right: B° — > ji~ . 
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Table 2: Expected combinatorial background events, expected peaking (B^ — > h + fi~) 
background events, expected signal events assuming the SM branching fraction prediction, 
and observed events in the B° s — > search window. 
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upper limits for the branching fractions are evaluated to be 

B{B Q s ^n + n~) < 1.3(1.6) x KT 8 at 90% (95%) CL, 
6(B°-^/iV) < 3.0(3.6) x 1(T 9 at 90% (95%) CL. 
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The B(Bg ) and B(B° — > n + n ) upper limits have been combined with those 

published previously by LHCb [6] and the results are 

B(B°-> / uV~)( 2 010 + 2011) < 1.2(1.4) x 10~ 8 at 90 % (95 %) CL, 
£(£°->-// + /O(2010 + 2011) < 2.6(3.2) x 10~ 9 at 90% (95%) CL. 

The above 90% (95%) CL upper limits are still about 3.8 (4.4) times the SM branching 
fractions for the B° s and 26 (32) times for the B°. These results represent the best upper 
limits to date. 
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Table 3: Expected combinatorial background events, expected peaking (B^ — > h + h' ) 
background events, expected B° — y signal events assuming the SM branching 

fraction, expected cross-feed events from B° s — y assuming the SM branching fraction 

and observed events in the B° — y [x + [x~ search window. 
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Table 4: Expected and observed limits on the B Q S — > fi + fi~ branching fraction for the 2011 
data and for the combination of 2010 and 2011 data. The expected limits are computed 
allowing the presence of B° s — > events according to the SM branching fraction. 







at 90% CL 


at 95% CL 


CL b 


2011 


expected limit 
observed limit 


1.1 x 10~ 8 
1.3 x 10~ 8 


1.4 x 10~ 8 
1.6 x 10~ 8 


0.95 


2010+2011 


expected limit 
observed limit 


1.0 x 10~ 8 
1.2 x 10~ 8 


1.3 x 10~ 8 

1.4 x HT 8 


0.93 



Table 5: Expected and observed limits on the B° — > branching fraction for 2011 

data and for the combination of 2010 and 2011 data. The expected limits are computed 
in the background only hypothesis. 







at 90% CL 


at 95% CL 


CL b 


2011 


expected limit 
observed limit 


2.5 x 10~ 9 
3.0 x 10~ 9 


3.2 x 10~ 9 
3.6 x HT 9 


0.68 


2010+2011 


expected limit 
observed limit 


2.4 x 10~ 9 
2.6 x 10~ 9 


3.0 x 10~ 9 
3.2 x 10~ 9 


0.61 
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